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Abstract 

Observations of topological defects associated with Stone- Wales-type transformations (z.e., bond 
rotations) in high resolution transmission electron microscopy (HRTEM) images of carbon nanos- 
tructures are at odds with the equilibrium thermodynamics of these systems. Here, by combining 
aberration-corrected HRTEM experiments and atomistic simulations, we show that such defects 
can be formed by single electron impacts, and remarkably, at electron energies below the threshold 
for atomic displacements. We further study the mechanisms of irradiation-driven bond rotations, 
and explain why electron irradiation at moderate electron energies (^100 keV) tends to amorphize 
rather than perforate graphene. We also show via simulations that Stone- Wales defects can ap- 
pear in curved graphitic structures due to incomplete recombination of irradiation-induced Frenkel 
defects, similar to formation of Wigner-type defects in silicon. 
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INTRODUCTION 

Stone- Wales defect [l, 2] - 5^1^(55-77) - is the simplest example of topological disorder 
in graphene and other sp^-hybridized carbon systems. It can be formed by rotating a C-C 
bond by 90° with regard to the midpoint of the bond - referred to as the SW transformation 
- so that four hexagons are turned into two pentagons and two heptagons. This defect 
has received considerable amount of attention , because it has the lowest formation 
energy among all intrinsic defects in grapheme systems, and because it presumably plays 
an important role in plastic deformation of carbon nanotubes (CNT) under tension [9] by 
relieving strain jlo|. It can also act as a source for dislocation dipoles 0, 11]. 



Regardless of being the lowest energy defect in graphene [12I and other sp^-hybridized 
carbon nanostructures, the 51^(55-77) needs about 5 eV to appear in graphene ^sl, IgI, 
and 3-5 eV in CNTs with a diameter above 1 nm 0, 15I , which should lead to a negligible 
equilibrium concentration of such defects at room temperature. However, recent advances 
in HRTEM have allowed the identification of such defects in graphene 13, 14] and CNTs js]. 
Moreover, SW transformations play an important role in the response of graphene to electron 
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irradiation p^, 115| , leading to changes in the morphology of vacancy- type defects [16[] and to 
their migration. Such changes are equally surprising, because the barrier for bond rotation is 
about 5 eV 0] , which should exclude thermal activation as a cause for S W transformation 
at room temperature during experimentally relevant time scales. Regarding irradiation 
effects, previous simulations [3] showed that an energy of ^ 30 eV must be transferred 
to a C atom in graphene in the in-plane direction for a bond rotation to occur. Also this 
cannot explain the frequently observed SW transformations under the usual TEM imaging 
conditions, since with typical acceleration voltages (< 300 kV) the transferred kinetic energy 
in the direction almost perpendicular to the electron beam will remain significantly below 
10 eV. 

Here, by combining aberration-corrected (AC-) HRTEM with atomistic computer simula- 
tions, we show that topological defects associated with the SW transformation can be formed 
in 5]9^-hybridized carbon nanostructures by impacts of individual electrons at energies even 
below the threshold for a carbon atom displacement. We further study in detail the mech- 
anisms of irradiation-driven bond rotations for pre-existing vacancy-type defect structures 
and how they transform and migrate due to electron impacts. At the same time we explain 
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why electron irradiation at moderate energies 100 keV) tends to rather amorphize [iJ] 
than perforate graphene. We also show via simulations that the 5'VI^(55-77) can appear 
in curved graphitic structures due to "incomplete" recombination of irradiation-induced 
Frenkel defects, reminiscent of the formation of Wigner-type defects in silicon jl9]. 



METHODS 



Experimental Methods 

Graphene membranes used in our experiments were prepared by mechanical exfoliation of 
graphite on Si/Si02 substrates and transfer of the resulting graphene sheets onto TEM grids 
as described previously {2^. For TEM experiments we used an FEI Titan 80 — 300 equipped 
with an image-side aberration corrector, operated at 80 kV. The spherical aberration was 
set to 15 fim and images were recorded at Scherzer defocus. The extraction voltage of the 
source was reduced to 2 kV and the condensor lens C2 was switched off in order to minimize 
the energy spread. Under these conditions, dark contrast in the images can be directly 
interpreted in terms of the atomic structure. Image sequences were recorded on a CCD 
camera with exposure times of 1 s and intervals of approximately 2 s. 



Computational Methods 

We carried out atomistic computer simulations based on the spin-polarized density func 



tional theory (DFT) implemented in the plane wave basis set VASP code [2l|]. The projector 
augmented wave potentials [22^ were used to describe the core electrons, and the generalized 
gradient approximation of Perdew, Burke and Ernzernhof [23] for exchange and correla- 
tion. We included plane waves up to a kinetic energy of 300 eV. The k-point sampling of 



the Brillouin zone was performed using the scheme of Monkhorst-Pack [241 the peri- 
odic dimensions. Structure relaxation calculations were combined with molecular dynamics 
(DFT-MD) simulations with a lower kinetic energy threshold and fewer k-points. 

Due to the high computational cost of the DFT-MD method, only a few simulations were 
carried out at this level. Whenever statistics needed to be gathered, we calculated the forces 



using the non-orthogonal DFT-based tight binding (DFTB) method [25|]. The main results 



were checked against DFT-MD. In total, we carried out ^ 27.000 dynamical DFTB-MD 



simulations. The simulated structures consisted of 120^200 atoms and were fully optimized. 
For the displacement threshold simulations, one of the atoms was assigned a kinetic energy 
T with the initial velocity vector pointing to a preselected direction. The initial temperature 
of the system was set to 5 K, although we observed no differences when carrying out the 
simulations for initially stationary atoms. Displacement threshold (minimum kinetic 
energy required to eject t: 
earlier DFTB results |2£ 



le atom) was found to be 22.50 eV, in a good agreement with 
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27|. It is also close to the DFT value (22.03 eV) [28]. For the 



annihilation simulations, various system temperatures were studied (500-1200 K) both to 
fasten the migration of the adatoms and to understand the effect of an elevated temperature 
(as will be mentioned below). 



RESULTS AND DISCUSSION 



Stone- Wales Defects due to Single Electron Impacts 

We begin the presentation of our results with the description of experimental observations 
of S'VI^(55-77) in HRTEM images. Several long image sequences, typically containing 
hundreds of images from clean and initially defect-free graphene membranes, were recorded 
at 80 kV. Occasionally, S'VI^(55-77) defects appear in individual exposures, as in the example 
shown in Figure [lb (Figure [lb with structure overlay). Remarkably, in most of the observed 
cases, isolated aSVI^(55-77) appeared in pristine graphene for one 1 s exposure, only to 
disappear in the following frame. Hence, the lifetime of this defect under the 80 kV electron 
beam in terms of irradiation dose is of the order of 10'' e~/nm^, the dose used for a single 
exposure. 

To understand the appearance and disappearance of S'VI^(55-77) , we carried out atomistic 
simulations for individual displacement events under the electron beam. After calculating 
T(i, we extended this calculation to all in-plane (0) and out-of-plane angles in the range 
9 G [0°,25°], Figure [21 Displacements with 6 > 25° would result in transferred kinetic 
energies of more than 2 eV below for an electron beam even at 120 keV. Since the 
displacement threshold increases for increasing ^, it is unlikely that this restriction would 
lead us to miss any significant electron beam-induced structural changes, especially for 
electron energies similar to those used in this study (80 keV). The calculated displacement 
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Figure 1: HRTEM images of an initially pristine graphene sample (a), followed by a 51^(55-77) 
(b). Frame (c) shows the same image with a structure overlay. 51^(55-77) disappears in the 
following frame (d). Scale bar is 1 nm. (See also video SI in Ref. [29].) 
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Figure 2: Relative displacement threshold T^'^ /T^^^ as a function of the displacement space angle 
(6^, (/)). Crosses mark the angles for which we observed impact-induced SW transformations. 

thresholds are shown in Figure [2] in a relative scale along with the space angles for which 
we observed the formation of aSVI^(55-77) . 

It is evident from Figure [2] that the SW transformation is a very likely event at dis- 
placement angles slightly away from the graphene plane normal {9 > 2.5°). The transferred 
kinetic energies (T) required for this process are below the displacement threshold for the 
corresponding cj) and 9 {Tp^) since no actual removal of the recoil atom is required for the 



5 



Figure 3: Typical SW transformation processes for low out-of-plane displacement angles {9 < 25°): 
(a) The "circle" process, and (b) the "nudge" process. The black spheres correspond to the recoil 
(displaced) atoms, and the mono-vacancy structure is highlighted with numbered carbon rings in 
the second panels. The last panel shows schematically the route of the displaced atom. (See also 
videos S8 and S9 in Ref. 

bond rotation to occur. Typically, T ^ T^^ — 1 eV resulted in the aSVI^(55-77) forma- 
tion, although for some space angles even T ^ T^'^ — 2 eV was enough. The probability 
for 5^1^(55-77) formation is particularly high for certain space angles, which is related to 
different mechanisms of SW transformation, as described below. 

The above- presented result is in clear contrast with the earlier simulation results for 

n 

graphite where no 51^(55-77) formation was observed for low 9. This discrepancy is 
caused by the neighboring graphene planes in the case of graphite: The displaced atom gets 
attached to the adjacent layer and does not therefore initiate a bond rotation. In Figure [3] we 
show the two processes which account for the majority of the SW transformations observed 
during our simulations. In the "circle" process (Figure[3K), the displaced atom circles around 
its neighbor, whereas in the "nudge" process (Figure [3h) it nudges the neighbor to cause the 
bond rotation. Note that the example cases are for the same (j) and almost same T, but for 
different 9. The resulting process for each displacement is an interplay of all three variables 

Similar mechanisms also exist for CNTs. However, as two new parameters (tube diameter 
and chirality) should be introduced for quantitative analysis of SW transformation, we did 
not study this process in nanotubes at length due to unreasonably high computational cost. 
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Figure 4: Experimental image sequence of a vacancy annihilation. Panel (a) shows a Vi(5-9) mono- 
vacancy, with an overlay in the inset, and (b) shows the same region in a later exposure, with no 
defect visible. Scale bar is 1 nm. (See also video S7 in Ref. [^].) 

Annihilation of Vacancies and Adatoms with Possible Formation of 51^(55-77) 

Since we frequently observed formation of vacancy-adatom pairs (adatoms play the role of 
interstitials in graphene and CNTs) in our simulations of electron impacts onto graphene and 
earlier in CNTs j26|, we also explored another possible mechanism of S'VI^(55-77) formation, 
which is based on "incomplete" annihilation of a Frenkel defect. This study was motivated 
by the peculiarities of the recombination of such a defect in bulk silicon. In that covalently- 
bonded material the recombination can give rise to either annihilation of the defect and 
restoration of the perfect crystal lattice, or to a Wigner-type defect Such topological 
defects are imperfections in the crystal lattice with the locally "correct" number of atoms 
(as opposed to vacancies and interstitials), with the atomic configuration separated from 
the perfect structure by a finite potential barrier. Such defects are deemed to also exist in 
graphite and aSVI^(55-77) can clearly be classified into this group. 



While carbon adatoms on graphene |32l, l33|] and CNTs |3J, l35| (especially those inside 



nanotubes) are mobile at room temperature, they can easily find vacancies in the system 
and annihilate. Indeed, we occasionally observed disappearance of vacancies in HRTEM 
image sequences. Figure [4] shows an example of a mono- vacancy that disappears during 
observation. This proves that mobile carbon atoms are present under our experimental 
conditions, and may recombine with vacancy- type defects. However, we never noticed the 
creation of a 5^1^(55-77) after an observed mono-vacancy. 
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Figure 5: (a) 51^(55-77) creation process as an outcome of a vacancy-adatom annihilation in a 
(6,6) armchair nanotube. (b) Carbon dimer sputtering process in a (10,0) zigzag nanotube. The 
black spheres stand for adatoms and the gray ones denote the dimer to be sputtered [marked in 
the first panel of (b)]. The nanotube structures are shown from inside the tube. The tube axis is 
in the horizontal direction. (See also videos SIO and Sll in Ref. [29.].) 

To simulate the annihilation process, we created nearby Frenkel defects (separated by 
a few Angstroms) in a graphene layer and small (6,6) zigzag and (10,0) armchair CNTs 
(radii r 4.1 A and ^ 3.9 A, respectively). We then heated the structures and collected 
statistics on the evolution of each system by running dynamical atomistic simulations at 
various temperatures (500-1200 K). Three possible outcomes emerged from the simulations: 
(1) perfect annihilation to the pristine structure (similar to the experimental images in 
Figure Sj), (2) formation of a S'VI^(55-77) (Figure [5k) and, surprisingly, (3) sputtering of a 
C2 dimer with a remaining reconstructed di- vacancy ^2 (5-8-5) (Figure [5b). 

For graphene, we always observed perfect annihilation in accordance with the experi- 
ments. However, if the adatom in graphene was placed on top of one of those two bonds 
in the nine-membered carbon ring which are right next to the pentagon, a aSVI^(55-77) 
was spontaneously formed without an energy barrier. Thus aSVI^(55-77) may also form in 
graphene due to recombination of Frenkel defects, but the probability for this process must 
be much lower than in CNTs. 

For (10,0) CNT, we obtained perfect annihilation in approximately 54% of cases, 5'VI^(55- 
77) was formed in approximately 34% of cases and dimer sputtering occurred in approx- 
imately 12% of cases. For (6,6) CNT, the values were 53%, 42% and 4%, respectively. 
With increasing temperature the probability to sputter a dimer showed a slight tendency 
to increase. We also ran the calculations for a (8,8) CNT (r ^ 5.4 A) at 800 K in order to 



estimate the curvature dependency of the results. The values did not significantly diflFer from 
those for the (6,6) CNT, except for a somewhat increased tendency to perfect annihilation 
and decreased sputtering (with probabilities of 56% and 2%, respectively). 

The formation energy of a spatially separated Frenkel defect in graphene within the 
DFTB model is approximately 11.1 eV. Transformation to form SW defect from this initial 
setup leads to an energy gain of 5.4 eV. By sputtering a dimer, graphene would instead gain 
0.4 eV. Hence, all three observed outcomes are energetically reasonable also for graphene. 
For nanotubes, the formation energies of both S'VI^(55-77) and a di- vacancy are lowered due 
to the curvature and stronger C-C bonds at pentagons [s^l. The corresponding energy gains 
are also higher, which can explain why the probability for aSVI^(55-77) defect formation 
and dimer evolution is higher in curved carbon nanostructures. The actual energies depend 



on the local curvature. It is also plausible that t 
in shrinking fuUerenes under electron irradiation 



le dimer evaporation process plays a role 



Stone- Wales Transformations in Vacancy- Type Defects 

S'VI^(55-77) defect represents the elementary case of a topological change in the graphene 
structure, i.e., a single bond rotation in the otherwise perfect structure. More abundant, 
however, are changes in the atomic configuration through bond rotations in the reconstructed 
vacancy defects, as recent experiments indicate Q,^. In presence of a (multi-) vacancy, the 
atomic configuration of a defect can be transformed between different metastable structures 
via bond rotations. It was presumed [14] that such SW transformations of vacancy-type 
defects would be stimulated by electron impacts, but the actual atomistic mechanism has 
not been hitherto unraveled. 

In order to get microscopic insight into irradiation-stimulated bond rotations near 
vacancy-type defects in graphene, we carried out a set of experiments and dedicated sim- 
ulations aimed at assessing the probability of SW transformations in the defect structures. 
In our experiments, we initially generated "defective" graphene by brief 150 kV electron 
irradiation, and then recorded image sequences of di-vacancies using 80 kV AC-HRTEM. As 
observed previously jl4|, the vacancies can transform between diflFerent configurations under 
the infiuence of the 80 kV electron beam. Moreover, the di-vacancies migrate and transition 
between the different reconstructed configurations via SW transformations, typically until 
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Figure 6: Experimental images illustrating SW transformations in the atomic structure of di- 
vacancies and migration of these defects. Panels (a-d) show sequential HRTEM images of the same 
defect, in the lower row with structure overlay, (a) (5555-6-7777) transforming into (555-777) 
(b), and V2 (5-8-5) (c). Each of these transitions can be explained by a single bond rotation. In a 
later frame (d), the defect is again a 1/2 (5555-6-7777), but shifted by one lattice parameter. Scale 
bar is 1 nm. (See also videos S2-4 in Ref. [29].) 

I — I 

they cluster into larger defects (see videos S5 and S6 in Ref. [29]). 

In Figure [HI we show an example of a di- vacancy defect that transforms between different 
reconstructed shapes [1^2 (5-8-5), 1/2(555-777), 1^2 (5555-6-7777)] under the electron beam. The 
changes in the atomic structure of these di- vacancy configurations can be described by SW 
transformations at the defect. Moreover, multiple transformations allow migration of the 
di-vacancy. Similar to SW{55-77) defect formation discussed above, the activation energy 
for these transitions is far too high to allow a thermally activated process with an observable 
rate at room temperature. Hence, on the basis of the observations, the activation energy for 
the transition must be provided by the electron beam. 

In order to confirm that the SW transformations at these defects are caused by single 
electron impacts, we carried out atomistic simulations of such impacts onto atoms near di- 
vacancies. Due to the computational cost related to many non-equivalent atoms present 
in the system and a large number of possible atomic configurations, we could not repeat 
the detailed analysis of the role of initial space angle of the displacement similar to pristine 
graphene, and therefore limited our simulations to the ^ = 0° case for all non-equivalent 
atoms at reconstructed di-vacancy structures. However, because the defects break the sym- 



10 



Figure 7: Changes in the atomic structures of reconstructed di-vacancy defects through SW trans- 
formations upon atom displacements perpendicular to the graphene sheet {0 = 0°); (a) V2(5-8- 
5)^ 1/2(555-777) and (b) V2 (555-777)^ ^2 (5555-6-7777). The black spheres indicate the recoil 
(displaced) atoms. Structurally equivalent atoms to the displaced ones are marked with circles in 
first panels. In the last panels, the circled atoms are those which could cause a backward transfor- 
mation upon displacement in addition to the recoil atom. (See also videos S12-14 in Ref. Q.) 

metry of the lattice, a directional preference for the displacements arises (atom displaced 
in the perpendicular direction will change its direction due to local strain). This effect is 
strong enough to facilitate bond rotations in reconstructed di-vacancy structures. In Fig- 
ure [71 we present examples of such processes for both V2(5-8-5)^ 14(555-777) and V2(555- 
777)^ V2 (5555-6-7777) transformations. 

Within these simulations, we never observed 1/2(555-777)^ 1/2(5-8-5) transformations. 
Because the symmetry of the 1/2(555-777) defect around the middle atom is the same as 
that of pristine graphene (the middle atom is represented as a black sphere in the last 
panel of Figure Uh) , one would also expect that a 7^ 0° displacement is required for this 
transformation, similar to the SW{55-77) case, although in principle the surrounding 
atoms could also cause this transformation. We noticed that the most likely di-vacancy 
transformation, at least for ^ = 0°, is the 1/2(5555-6-7777)^ 1/2(555-777). The 1/2(5-8- 
5)^ 1/2(555-777) was the least likely one of those observed. We never observed a 1/2(5-8-5)^ 
2x(5-7) transformation fl^ during our simulations, which we also attribute to the limited 
simulated conditions {9 = 0°). Curiously, however, we did observe one transformation in 
which a 1/2(5-8-5) di-vacancy directly migrated one step in the zigzag lattice direction. 
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Another interesting observation originating from these simulations is the fact that the 
displacement threshold for atoms in the central part of the reconstructed defects [V2(555- 
777) and V2 (5555-6-7777)] are higher than that for pristine graphene (by as much as 5%). 
This may explain why defect structures tend to grow into larger and larger amorphous 
patches instead of collapsing into holes under continuous electron irradiation at low voltages 
(< 100 keV) fl^: Even when atoms are removed from the defected area, the displacements 
occur at the edges of the existing defects rather than at the central part where the local 
atomic density is already lower. Clearly, since the core structure of these defects consists of 
carbon hexagons, there must exist a limiting size above which the displacement threshold 
becomes similar to that of ideal graphene. 



CONCLUSIONS 



To conclude, by combining AC-HRTEM experiments and atomistic simulations, we have 
shown that the bond rotations which lead to creation of topological defects in carbon nanos- 
tructures are caused by single electron impacts or incomplete annihilation of Frenkel defects. 
This explains the discrepancy between experimental observations of Stone- Wales defects and 
their relatively high formation energy and even higher energy barrier for bond rotation. 

The SW transformation in graphene can be initiated at least in two different ways upon 
electron impact (involving a "circling" or "nudging" motion), and for almost any space an- 
gle, provided that enough energy is transferred from the electron to the target atom. Our 
simulations indicate that S'VI^(55-77) can appear as a result of "incomplete" recombina- 
tion of a Frenkel defect reminiscent of the formation of Wigner-type defects in silicon |l9 |. 
However, this is much more likely in the case of local curvature, as in nanotubes. More 
surprisingly, we also observed sputtering of C2 dimer as a result of annihilation of a Frenkel 
defect in carbon structures with high curvature. 

Moreover, we showed that the displacement threshold of atoms in the central area of 
reconstructed defects is higher than that of pristine graphene, by as much as 5%, which 
explains why defected graphene under low-energy electron irradiation (< 100 keV) tends 
to turn graphene into a two-dimensional amorphous structure jl4| instead of a perforated 
membrane. For different di- vacancy structures, even displacements in direction perpendicu- 
lar to the graphene layer can initiate SW transformation and thus local structural changes 
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and defect migration. 

Our results provide microscopic insight into the irradiation-induced changes in the atomic 
structure of carbon nanosystems under electron irradiation, and taking into account the 



37, 



interesting electronic properties of defects associated with SW transformations [IJ, 

n 

may open new avenues for irradiation-mediated engineering [39^ of carbon nanostructures 
with next-generation electron microscopes. 
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